Kondo Resonance of a Microwave Photon 



Karyn Le Hur^'^ 

Centre de Physique Theorique, Ecole Polytechnique, CNRS, 91128 Palaiseau Cedex, France 
^Department of Physics, Yale University, New Haven, CT 06520, USA 

We emulate renormalization group models, such as the Spin-Boson Hamiltonian or the anisotropic 
Kondo model, from a quantum optics perspective by considering a superconducting device. The 
infra-red confinement involves photon excitations of two tunable transmission lines entangled to an 
artificial spin-1/2 particle or double-island charge qubit. Focusing on the propagation of microwave 
light, in the underdamped regime of the Spin-Boson model, we identify a many-body resonance 
where a photon is absorbed at the renormalized qubit frequency and reemitted forward in an elastic 
manner. We also show that asymptotic freedom of microwave light is reached by increasing the input 
signal amplitude at low temperatures which allows the disappearance of the transmission peak. 
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The asymptotic confinement phenomenon in the infra- 
red limit is omnipresent in condensed-matter systems and 
it plays a crucial role in quantum impurity systems, such 
as the Kondo model describing a single spin-1/2 parti- 
cle interacting with a bath of conduction electrons [1]. 
The Kondo effect can also be considered as an exam- 
ple of asymptotic freedom, i.e., the coupling of electrons 
and spin becomes weak at high temperatures or high 
energies. This model introduced to describe resistance 
anomalies in metals with magnetic impurities embodies 
the "hydrogen atom" of many-body physics [2, 3]. Dis- 
tinct aspects of this infra-red confinement phenomenon 
can also be addressed through a one-dimensional boson 
bath (transmission line) entangling a spin-1/2 particle or 
two-level system resulting in the Spin-Boson model which 
can be mapped onto the anisotropic Kondo model and 
exhibits a plethora of interesting phenomena such as an 
underdamped-overdamped crossover in the spin dynam- 
ics and a quantum phase transition [4-6]. In this Letter, 
we consider the superconducting Josephson circuit of Fig. 
1, which allows to investigate the quantum entanglement 
in the Spin-Boson model and therefore properties of the 
anisotropic Kondo model through transport of photons. 
In the underdamped limit, we prospect to reveal a many- 
body Kondo resonance in the elastic power of a transmit- 
ted microwave photon. This circuit offers the opportu- 
nity to export many-body physics in quantum optics. 

The superconducting system comprises an artifi- 
cial spin or double-island charge qubit [7-9] interact- 
ing with the zero-point fluctuations of two long one- 
dimensional transmission lines envisioned from tunable 
one-dimensional Josephson junction arrays [10-12]. In 
order to maximize the elastic transmission of a microwave 
photon, the spin-1/2 object is built from a superconduct- 
ing double Cooper-pair box where spin up and spin down 
states refer to the two degenerate charge states (0, 1) 
and (1,0), respectively corresponding to one additional 
Cooper pair on either island [13]. Recently, geometries 
involving artificial atoms and transmission lines or cavi- 
ties have already been realized experimentally [14, 15]. 
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FIG. 1: (Color online) Superconducting circuit envisioned 
from a double-island charge qubit coupled to two one- 
dimensional Josephson junction arrays allowing to produce 
a Kondo resonance in the elastic power of microwave light. 



Below, we will assume that the charging energy corre- 
sponds to the most dominant term in the Hamiltonian. 
In fact, close to a charge degeneracy line [7], we can apply 
the pseudospin representation for the charge states (0,1) 
and (1,0) reinterpreting them as spin- up and spin-down 
eigenstates of the operator Uz [13]. The effective detun- 
ing e = {EiQ — Eqi) — )• 0, where Eio (E'oi) corresponds 
to the energy of the spin-down (spin-up) eigenstate, can 
be adjusted through the gate voltages Vga and Vgb. 

Transfer of Cooper pairs between superconducting is- 
lands and leads is described through the Josephson terms 
EjL and Ejfi in Fig. 1. In the weak tunneling limit 
{Ejl.Ejr) < min(£:ii - Eio^Eoq ~ Eia) one can per- 
form a standard perturbation theory and cotunneling of 
Cooper pairs then costs an effective energy [13]: 
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where En {Eqo) corresponds to the energy to add 
(remove) one extra Cooper on the double-island. 
The Josephson Hamiltonian then takes the form 
-{Ej/2)a+ exp[i{(f>i - (f>r){x = 0)] + h.c. [16] where the 
Josephson phases (j>i {x ~ 0) and (pr {x = 0) of the left and 



right one-dimensional transmission lines read {j = I, r 
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This Joscphson term captures the transmission of a given 
Cooper pair across the system in Fig. 1. The su- 
perconducting reservoirs are explicitly modeled by one- 
dimensional transmission lines revealing low-energy pho- 
ton excitations. The left and right transmission lines are 
described by two distinct sets of harmonic oscillator (pho- 
ton) operators bik and brk- Below, we consider the limit 
where (CL,Ci?) < and {Ejl,Ejr) < H^/Lt in Fig. 
1. The spatial solution of the modes can be expressed 
in terms of the wavevectors k w rmr/{2£), where m is 
odd for symmetric modes and even for antisymmetric 
modes, and a tra.nsinission line is diagonalizcd introduc- 
ing bosonic creation and annihilation operators. Here, £ 
corresponds to the length of each transmission line and 
c = C(/a to the capacitance per imit length; a is the size 
of a unit cell in each transmission line and we consider 
the thermodynamic limit a/C — >■ 0. The photon waves 
propagate at the speed v = uj^a where = ^/^/LfCt, 
the inductances are defined in Fig. 1, and = v\k\. 

To build an explicit analogy with the spin-boson 
Hamiltonian, we rewrite the Joscphson term as a trans- 
verse field Hj = —[Ej/2)ax performing a unitary trans- 
formation or spin rotation (see footnote in [29]). Such a 
procedure, also referred to as a polaron transformation 
[4, 5], has been applied in the case of a spin- 1/2 interact- 
ing with the sound modes of a Bose-Einstein condensate 
[17]. Since the Hamiltonian of the transmission lines does 
not commute with the spin rotation this produces an ef- 
fective interaction between the two-level system and the 
photon excitations. This term can be combined with the 
capacitive couplings Cl and Cr of Fig. 1. More pre- 
cisely, since the electrical potential (operator) at the end 
of a transmission line, i.e., at x = 0, takes the form 

Vj = 0) = ^ E Vjy^kibjk + b]k), (3) 
'^'^^ k>0 

this results in the Spin-Boson Hamiltonian: 
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Following Ref. 7 and the notations of Fig. 1, wc have de- 
fined the total capacitances seen by each superconducting 
island: C^a = CL + Cga + Cm and Ce^ = CR + Cgb + Cm- 
The analogy with the Spin-Boson model [4, 5] becomes 
complete when rewriting the Hamiltonian in terms of the 
symmetric and antisymmetric bosonic combinations: 



bsk = cos6'6;fc -I- sin6'6rfc 
bak = sin ebik - COS 6brk- 



(6) 



Choosing cos = jr/ \/ if + Ir sin = ^i/ \/7r + if, 
we note that the boson operator bak only couples to the 
two-level system through the coupling = ak\Ji1 ^" 7r ■ 
Each transmission line mimics a physical resistor then 
producing dissipation in the system. In the present cir- 
cuit, the spectral fimction of the environment is defined 
as J(w) = {'KlK)Y,k>Q^kK^ - ^k) = 27rfiau;e""/"<^ 
where Wc ^ Ejjh. represents the high-frequency cutoff 
of this Ohmic environment [30] and a is given by 
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Here, Rq = h/{2e)'^ denotes the quantum of resistance 
and R = ^jLtjCt is the resistance of each transmission 
line. It is instructive to observe that in the limit of negli- 
gible capacitances Cl and C/f the system naturally con- 
verges towards the symmetric condition 7; = 7^ = — 1. 

The Spin-Boson Hamiltonian with Ohmic dissipation is 
intimately related to the Kondo model in the anisotropic 
regime via bosonization [18]. Other Spin-Boson Hamilto- 
nians such as the Jaynes-Cummings model, in contrast, 
involve a two-level system interacting with a single mode 
of a cavity [19]. Other impurity models with photons 
have also been considered [20, 21]. We are interested in 
the underdamped regime (0.1 < a < 0.2) of the Spin- 
Boson model where the two-level system displays visible 
Rabi oscillations but dissipation modifies the qubit fre- 
quency which is related to the Kondo energy [4, 5] 
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To understand the physical content of the energy Er, 
it is relevant to apply the unitary transformation U = 
exp(^/ - Ar) where Aj = J2k>o l^(^lfe - bjk)crz/2 such 
that the Hamiltonian can be rewritten as {H = WHU): 



The charge operators on the two islands take the forms H 
Q(, = ^(1 o-z) and Qa = ^(1 - (Jz)- Hereafter the 
detuning will be fixed to e — >■ and ak = (2e/-\/ cC)\/hu}k. 
The couplings 7^ and 7; are given by: 
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where the phases = —ji(f)i{x = 0) and $r = 
—"Iriprix = 0) contain Joscphson physics as well as 
(weak) charging effects. Then, we can define an effec- 
tive transverse field acting on the dissipative two-level 
system as A = Ej{cos{^i — $r)) such that the ar- 
tificial atom is described by the effective Hamiltonian 
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Heff = — (e/2)(T^ — (A/2)(Ta;. Bethe ansatz calculations 
[22, 23] and the adiabatic renormalization [4] in the un- 
derdamped Hmit where 0.1 ^ a ^ 0.2 indeed confirm 
that A = Efj. The bare qubit frequency Ej/H of the 
two-level system is modified due to the strong renormal- 
ization effects associated with the photon bath. One way 
to experimentaUy measure the Kondo energy i^jj would 
be through charge measurements since the Fermi liquid 
ground state imposes that (a^) cx e/En at small detun- 
ing and low temperatures fc^T <C Er, and the prefactor 
is accessible from Bethe Ansatz calculations [23]. 

Below, we show that in the underdamped regime and 
for temperatures ksT <^ Er, the Kondo energy Er 
can be directly measured based on the (elastic) resonant 
propagation of a photon. When the system is driven by 
an external coherent source, the drive, the circuit and the 
outgoing waves can be treated through the input-output 
theory [24] . Previous works have studied the limit a 
where many-body effects can be fully ignored (the elas- 
tic resonance is centered at the bare frequency of the 
two-level system and converges to a ^-function) [14]. We 
assume perfect transmission of the microwave signal in 
the transmission lines such that the input signal reads 
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Here, to < t denotes a time in the distance past before 
any wave packet has reached the two-level system. Simi- 
larly, an output field in the left transmission line at time 
ti > t being a time in the distant future after the input 
field has reached the double-island Cooper box system 
reads 

Vr\t) = E it {e-'^'^'-'^'^bikiti) + e'^''^''''^bl{t,)^ . 

k<0 

(11) 

Through the Heisenberg relation bik = {i/h)[H,bik] — 
—iidkbik + {i/'2h)^iakO'z we relate the properties of the 
input signal to those of the two-level system. 

Below, since we focus on the underdamped limit 
of the Spin-Boson model which is characterized by a 
(Rabi) resonance at uj = lor, we establish (^^(aj)) w 
7;x(a;, Pin){V"^{ijj, Pin)) for frequencies in the vicinity of 
ujR where Pin = {{V'l^'^)'^) / R is the average input power 
[31]; see EPAPS [25]. Then, we can introduce the reflec- 
tion coefficient r(w,P,„) = (l/™*(a;, P„))/Ct^"(a;, P„)). 
Defining the output signal in the right transmission line 
as 



yout ^ ^ g (e-'"'=(*'*i)&,fc(ii) + e^"''(*-*^)6|;fc(<i)) . 

fc>0 

(12) 

the transmission coefficient is t(uj,Pin) = 
(y,°«*(a;,P,„))/(V^™(a;,P™)). Using the Heisenberg 
relation of bik with the definitions above we obtain (see 
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FIG. 2: (Color online) Normalized elastic transmitted power 
tt* (lj) as a function of frequency and driving power for 7; — 
7r. The parameters are chosen as a = 0.15, ljr = 1 — Pr, 
Ej ~ 1.9, ojc = 50 and h — 1 (we set the ratio Ej/fkOc to 
be moderate, but the "many-body" resonance frequency or 
renormalized qubit frequency ujr is distinct from Ej/h). 



EPAPS [25]): 
r(w. Pi, 
t{uj, Pi, 
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First, we consider the linear regime where the am- 
plitude of the input signal is very small. Pin ~^ (see 
footnote in [29]). For frequencies close to the confine- 
ment frequency, again assuming the underdamped limit 
(0.1 < a < 0.2), we derive an expression of the spin sus- 
ceptibility which agrees with Numerical Renormalization 
Group resuhs [26]; see EPAPS [25]. This leads to 



X(w) 



UJr/H 
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where the dissipation factor takes the form 7(0;) = 
ujfiJ{uj)/H and is in agreement with the (many-body) 
Fermi-liquid type ground state [5] . In the linear regime of 
small input power, we check that the scattering matrix is 
unitary, |r|^ -|- = 1, since J(a;fl;)9mx(ti-'i?,) = 1 show- 
ing that the photon propagation is purely elastic close 
to the resonance (see Fig. 2). We corroborate that the 
normalized (elastic) power tt*(ujR) flowing to the right 
transmission line reaches unity since here 7; « 7^. 

In the underdamped regime, the photon propagation 
across the system is characterized by a many-body res- 
onance at the frequency ujr [4-6] : a photon is absorbed 
at the frequency ujr and reemitted forward in a purely 
elastic manner. In the underdamped regime of the cir- 
cuit, the qubit is described by a resonance which turns 
the "photon-l-Cooper pair" system into an ideal conduc- 
tor. The phase associated with the reflection coefficient 
satisfies the following properties. For small 7/, the phase 
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vanishes since V°"^ = V"^ for an open termination and 
for 7^ = the phase is consistent with the Kondo-type 
2 X 7r/2 phase shift of a right-moving wave. 

In fact, the appearance of resonances in snch a circuit is 
not so surprising. For example, let us ignore the Coulomb 
blockade physics in the two islands completely and re- 
place the Josephson junctions Eji and Ejj^ by piirely 
linear inductances Ll and Lji. When Cl = Cr = C 
and Ll = Lr = L, then we corroborate a resonance 
with r = at the frequency w* = 1/^/CL + 2Cm.L; how- 
ever, we emphasize that ui* is distinct from which 
has a many-body origin. In addition, nonlinear effects 
unavoidably appear in the Josephson circuit of Fig. 1 
when increasing the amplitude of the input signal. Un- 
der a strong drive, this produces the accumulation of 
a macroscopic number of photons in the left transmis- 
sion line which will cause the saturation of the two-level 
system excitation and the destruction of the resonance 
peak. More precisely, evaluating the Franck-Condon fac- 
tor (cos(<l?i - = exp -[(($; - $r)2)/2] when increas- 
ing the input signal amplitude, we observe that {^f) 
yields an extra contribution (note, Pin = {Vi{x = 0)^)/R) 
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and the symbol ' refers to momenta such that utg ^ lur 
in the case of a monochromatic signal with a frequency 
close to u)r. The Josephson process in Eq. (9) becomes 
exponentially diminished and this results in an exponen- 
tial suppression of 3mx(w = cor) (see EPAPS [25]) 



'<smx{i^R, Pin)J{(^R) = exp - 
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where Pr = tvjj\. The nonlinearity of the two-level sys- 
tem produces an exponential decrease of the spin sus- 
ceptibility at cj = (jjR. Then, this causes the disappear- 
ance of the Rayleigh transmission resonance; see Fig. 2. 
When Pin ^ Pr one reaches the asymptotic freedom 
of microwave light where r{u)R) ~ 1: sec Eqs. (13) and 
(16). These nonlinear effects are driven by the Josephson- 
type Hamiltonian in Eq. (9). Increasing the driving 
power Pin, the scattering matrix becomes non-unitary 
since J{ujR)Qmx{^^R, Pin) < 1 (which hides the presence 
of additional inelastic Raman corrections [14]). An open 
question concerns the inelastic Raman spectrum, which 
should become prominent in the overdamped limit of the 
Spin-Boson model. The asymptotic freedom of light, re- 
sulting in r{u!R) ^ 1 can also be reached when increasing 
the temperature producing a prominent decoherence of 
the two-level system and a strong decrease of (a^) and of 
the Josephson coupling in Eq. (9); see EPAPS [25]. 

To summarize, in the underdamped regime of the Spin- 
Boson model, sending a microwave photon produces a 
many-body Kondo resonance. For moderate and acces- 
sible values of the dissipation strength 0.1 ^ a ^ 0.2, 



the confinement frequency lur is clearly distinguishable 
from the bare Josephson frequency Ej/Hoi the two- level 
atom as a result of the continuum of photon modes in 
the (very long) transmission lines. The width of the res- 
onance peak reflects the Fermi-liquid type ground state. 
We assumed that the detuning e and thermal effects 
through kgT are smaller than the Kondo energy Er, and 
in the underdamped regime, Er is not too small com- 
pared to the Josephson energy Ej. For Josephson junc- 
tion arrays with large resistances, this circuit would offer 
the opportunity to study overdamped spin dynamics and 
quantum phase transitions using nonlinear optics. Note 
that the development of techniques such as the numerical 
renormalization group [27] or a stochastic-type approach 
[28] would be necessary to combine the traditional input- 
output theory with many-body physics of quantum impu- 
rity models. Finally, such superconducting quantum de- 
vices can be used for controllable (single-)photon sources 
in which a plethora of novel effects related to many-body 
physics and nonlinear quantum optics can be realized. 
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We introduce the input-output theory and we derive an analytical expression of the dynamical 
spin susceptibility in the weak-coupling limit in the presence of an incoming signal which illustrates 
the concepts of the Renormalization Group and agrees with the Numerical Renormalization Group. 
We also investigate nonlinear effects close to the Kondo (resonance) frequency or renormalized qubit 
frequency ujr when increasing the amplitude of the input signal. 

PACS numbers: 03.65.Yz, 03.75.Lm, 42.50.-p, 85.25.-j 



INPUT-OUTPUT ANALYSIS 



First, wc apply the concepts issued from input-output theory for the left transmission line and we derive Eqs. (13) 
in the main text. Let us start with the Heisenberg equations of motion: 

bik = '^[H,bik] (1) 
i , 

Note that the original Hamiltonian and H defined in the main text are related through a unitary transformation 
analogous to U. Now, bik + bj/^ — {i/h)[UHU~^ ,bik + bj).] — {i/h)[H,bik + b]/.] with bik — dbik/dt where U has been 
defined in the main text. Therefore, bik — i'i/f'')[H,bik]- 

By integrating the Heisenberg equations of motion, this results in: 

bik{t) ^ eM~^^k{t - h)]bik{to) + ^liau I dre"^''^'-^^ a,{T) (2) 



blit) = exp[t^k{t-to)]bl{to)-^jiak /*dre-'=(*--)a,(r). 



(jjk is the bare frequency of the transmission line (when setting 7; = 0), and therefore we identify to < t as a time in 
the distance past before any wave packet has reached the qubit (two-level system) . Assuming perfect transmission of 
the signal in the left transmission line then we identify: 

y;"(i) = yz(2; = 0,i) = ^|*^ (expH^fe(i-io)]&/fc(io)+exp[ic<jfe(t~to)]6L(io)) • (3) 

A; 

Now, let us consider an output field in the left transmission line at time ti > t being a time in the distant future 
after the input field has reached the qubit. In this case: 

/ti-i-+oo ^ 
dr— 7zafce-''"'^(*-"V^(r) (4) 

ti->-+oo ^ 



b\^^{t) = expiwfc(t-ii)6|^(ii)+ / rfr— 7,afce'''"^(*-"V,(r). 



Then, identifying: 



Vr\t) = Vr{x - 0,i) = ^ g (eM-i^k{t - h)]bik{h)+eMi^kit~h)]b\^{t{)) , (5) 

k 



after Fourier transformation and comparing Eqs. (2) and (4), we get: 



(We check that the results are unchanged restricting the sum in the expression of VJ*" to X^g>o ^^'^ 
expression of VJ°"* to J2q<o- We can always change g — >■ — g > in the sum for VJ""*.) 
Defining {az{uj)) = x(w, Pin)"li{Vl"{u), Pin)) and combining the two lines then we find 

{Vr\i^,Pin)) = {Vr{iO,Pin)) (l + 2i^^J{i0)x{uJ,Pin)) ■ (7) 

V "I" Tr / 

The reflection coefficient in the main text then is defined as: 

(youts 

ri^^Pin) = ^ = l + 2^^^^^J{0J)x{uJ,Pi„). (8) 

Note also that crj(r) = azir) implies X*('^) = x{~<^)- 
Proceeding in the same way with brk then we obtain: 

t{0J,Pin) = -^^^J{w)x{0J,Pin)- (9) 



DYNAMICAL SUSCEPTIBILITY IN UNDERDAMPED AND LINEAR REGIME 

Here, we derive an expression of the dynamical spin susceptibility x(w) for frequencies co close to the confinement 
frequency lur based on a weak-coupling expansion assuming a negligible input power. We apply a procedure analogous 
to the Non-Interacting Blip Approximation which has been shown to be valid in the limit 0.1 ^ a ^ 0.2 [1-3]. 

Let us start with the Heisenberg equation of motion where the Hamiltonian H is given in Eq. (4) of the main text 
and evaluate (^^(cj)) {H and the original Hamiltonian are related through unitary transformation which commutes 
with the operator cr^): 

&z='^[H,az] = -^ay. (10) 



Then, we obtain: 



\ k>0 fe>0 / 



For simplicity, here we assume that e = strictly. From these two equations, we obtain: 



Ej . Ej I Ej 
az = —^ay = -—l — 



\ fe>0 fe>0 / / 



As a result: 



cTz + (^) = +1^ ( I E ^^(^l'^ + blk) - I E ^^('^rk + bl,)) a^. (13) 

^ ^ \ fe>0 k>0 J 



On the other hand, we get: 



(14) 



\"' k>0 k>0 I 

This is equivalent to: 

oAt) = oAh) - 1^ £ (I + - I E ^^^^rk + oziil). (15) 



fe>0 /c>0 
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During the integration procedure, ti < t represents an initial time which is set arbitrarily (but at t = we must 
satisfy that the incoming signal has already reached the two-level system since we want to study the spin response to 
an incoming microwave signal). This results in: 

^ ^ \ fe>0 fe>0 / 

\ k>0 k>0 J 



\ q>0 q>0 



. + bl){t')]a,{t'). 

q>0 '" q>0 



If y,'" — >■ 0, at the time ti, we can assume that the spin-boson model is in its ground state and therefore from Bethe 
Ansatz we rigorously identify [4] {axiU)) = (ax) — > Eji/Ej (again, assuming that we are in the weak-coupling regime 
with < a <C 1/3; for another derivation, see below). More precisely, we obtain the exact equation (for all a) [4, 5]: 

where 

_ e-V(^-^")r[l-l/(2-2a)] 
r is the incomplete gamma function and 6 = alna + (1 — a) ln(l — a). Then, we apply the mean- field decoupling: 

k>0 fe>0 



The mean-field decoupling (averaging) is usually valid in the underdamped limit of the spin-boson model [3] , and then 
the input signal V;*"(t) = VoCOs{Lot) then (simply) mimics a time-dependent detuning e(f) acting on the two-level 
system. This approximation is well justified for time scales smaller than the time ^ l/w/? which corresponds to the 
crossover scale toward the strong-coupling regime. Here, we extend this approximation for frequencies slightly below 
the confinement frequency. This will result in ^^(w) « 7;x('^)l^'"(w) in the vicinity of the confinement frequency. 
Now, let us focus on the (real part of the) last term in Eq. (16). Assuming ti — t^ — oo, this gives: 



it + 7r 



fe>0 

/* dt' E "fc coth cos(c^fe(i - t')){aS')), 



where (3 = l/fc^T. When Pin = ((V;™)^)/i? — Vq /2R — > 0, the boson modes are taken to be in thermal equilibrium. 
After Fourier transformation (assuming t ^ 0), then to compute the real part one needs to evaluate: 

and V denotes the principal part of the integral. 

In particular, for frequencies ui ^ Ej/h which corresponds to the resonance (Rabi) frequency of the non-dissipative 
spin and for T — ^ 0, this renormalizes the term [E j / hY {a zioj)) in w|j(cr^(u;)) where ujr = Er/H formally obeys: 

The present approach assumes a weak-coupling a for which w_r can be identified to the Kondo frequency: 

"K»^(#i)*. (23) 
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This formula is formally valid for kgT Ej since the dominant contribution stems from frequencies higher than 
Ej/h. The renormalization procedure in the Spin-Boson model, consisting to integrate out high-frequency photon 
modes, results in an effective Josephson energy Er < Ej [1]. 

To compute the imaginary part leading to dissipation (damping), first let us consider that ksT > i?/?, such that 
the weak-coupling decoupling in Eq. (19) becomes rigorous. For frequencies w ~ w^, then we find a damping term 

-n(a;,T)(a,(w)), (24) 

where 

7(u;,T) = M:j(^). (25) 

We check that 7 is odd in frequency which guarantees that X*(w) = x{~^)- Now, let us decrease progressively the 
temperature such that ksT Er. Then, the damping term takes the quantum form: 

7=^J(a;). (26) 

The last form is the correct form of the dissipation term in the quantum limit {ksT < En): This is reminiscent of 
the Korringa-Shiba relation at low frequency reflecting the Fermi-liquid Kondo ground state [6, 7]: 

^mx{uj) = J{uj)^ex{uJ = 0). (27) 

More precisely, for frequencies co in the vicinity of ujr, we find the following dynamical response: 

(a^iuj)) {-LJ^ « ^^ T^'^rH ^ (28) 

The spin dynamical susceptibility for u close to con then takes the underdamped form (see Eq. (14) in the main text): 

X(-) « , ^"f^^^. , (29) 

This form of xi^) is in fact in good agreement with Numerical Renormalization Group results for 0.1 ^ a ^ 0.2 [7] 
and therefore in the main text we will assume this form of the spin susceptibility for 0.1 ^ a ^ 0.2. Note that this 
results in the equalities: 

Qmx{ojR)JiLOR) = 1 (30) 
SRex(wfl) = 0. 

Combining these equations with Eqs. (13) in the main text, this shows that the scattering matrix is unitary when 
Pin — >■ for w ~ ujR. Most of the spectrum is elastic close to the resonance. 

Eq. (29) is the susceptibility used in the main text, for the underdamped regime (0.1 ^ a ^ 0.2) and not too low 
frequencies. Note that in the low-frequency limit w <C ujr, replacing J2k>o 

-akli{hk + b\k)(Tz/'2. by -V;'"7;£7^/2 would 

become less accurate. 



NONLINEAR EFFECTS WHEN INCREASING THE INPUT POWER 

Prom Eq. (19), Eq. (29) must be in fact understood as: 

, . Ej{a,{U))/h^ 

X(w) « 2 2 • f ^ ' (31) 
— — 

To investigate the effect of the input signal amplitude (power), we notice that the result {ax{ti)) — >■ Er/Ej can be 
easily understood in terms of the polaron transformation resulting in the transformed Hamiltonian H in the main text. 
The effect of the bath is included only through the modification cr+ — > cr+ exp i(<I); — $r) and cr_ U- exp — i((f>; — $r)- 
Therefore, at sufficiently weak couplings (a <^ 1/2) [4] and T = 0, this results in (ctx) = ((Tx)c«=o x (cos($; — $r)) = 
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(cos(<J>; — <&r))- Again, wc have assumed e 0. This is the exponential dressing Pranck-Condon factor [1]. This 
implies that Eq. (31) is formally equivalent to: 

Ej (cos($,-$,)(t,)) 
xM = ^^|_^2_i^(^) • (32) 

We identify (cos($; - $r)) = exp -[(($; - $r)^)/2]. Now, let us analyze the terms in (($; - ^r)^). One finds: 

When Pin and T — >■ 0, we get {bikb]^) = 1 for all k and similarly for {brkbl,.) (this is still true for finite a since 
(cTx) = 0), and using the adiabatic renormalization [1]: 

« + *?> = E7i^-?r S*. = 2„.„f|^V (34) 



Assuming a weak coupling a, this allows to recover the result from Bethe Ansatz [4] 

(cos($, - $,)) ^ En/Ej. (35) 

It is relevant to observe that when Pin and a <C 1/2, one can formally replace (cos(i>; — $r)) by 1 — (^^ -|- <^1)/2. 

Now, let us slightly increase the input signal amplitude; recall, V"/"(t) = Vbcos(a;t) and we are interested in the 
elastic propagation of a photon with frequency lo ^ ujr. The (time- averaged or mean) input power then takes the form 
Pin = Vq/2R = NhiV[j , where N = dN/dt and A'^ represents the number of photons at a; = in the left transmission 
line (with an energy ^ fii^n)- Note that following the notations of the main text, one can also identify: 

= E = Y.^.{vic){bib,,). (36) 

9>0 ^ ' 9>0 



p = 



Using the fact that for a monochromatic source 

y/" = ^ (exp(iwt) + cxp(-iwi)) (37) 

= E §('''«(*) + ^l(*))' 

and that for the incoming modes biq{t) = biq{0) exp{—iiOqt) then, in the sum above, this selects momenta q such 
that ojq ^ Lu ^ LOj^. Therefore, there is a novel contribution to {^f) stemming from ujq ^ ujr in addition to the 
high-frequency contribution in Eq. (34). Since the initial time ti is set arbitrarily (but, at this time the incoming 
signal has reached the two- level system), when performing different measurements one must also average over ti and 
therefore {^f) yields an extra contribution equal to 

Jf—y, E ".7, {b]qkq) = -^J—^, (38) 

and the symbol ' refers to momenta such that ujq ^ u ^ lur in the case of a monochromatic signal. This results in an 
exponential suppression of Ssmx(a; = ojr); see Eqs. (15) and (16) in the main text. Again, this expression assumes 
that frequencies of interest lie in the range of cjr. We have also assumed that photons with a (much) larger frequency 
are still in a thermal equilibrium. 

In the underdamped limit, one photon is perfectly transmitted in the time scale I/lor as shown in the main text. 
We thus define the associated power Pr = ujR{hjjR). The magnetic susceptibility in the frequency domain lo ljr 
then reads (assuming that P,„ < Ej/h; see below) 

/ D \ (wfi/;i)exp-^ 
X(W ~ OJR, Pin) w ^ , (39) 

where 

A=^-^n^f. (40) 

When increasing the input signal amplitude this will produce a macroscopic number of photons with an energy Klor: 
the saturation of the artificial atom excitation manifests itself in a substantial decrease of the photon transmission. 
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ASYMPTOTIC FREEDOM OF MICROWAVE LIGHT 



In Eqs. (21) and (34), when Pj„ > Ej/h, note that formally the low-frequency cutoff of the integrals must be 
changed into Pin/Ej {hPin/Ej becomes a large energy scale controlling the departure from {bikb\f^ = 1, i-e., from 
the ground state) and the characteristic frequency in Eq. (39) becomes 

= ^ (1 - aHEjujP^r,)) . (41) 

Finally, we discuss thermal effects. At finite temperatures, performing a thermal average, one finds that (a^) substan- 
tially decreases for fc^T ~ E^ [8], since the artificial atom lies in a highly mixed state, producing a strong suppression 
of the spin susceptibility (this reflects that the reflection coefficient of the microwave light is strongly enhanced). 
Further, for {j3hEj) <C 1, from Eq. (21) we also predict that the characteristic energy of the artificial atom takes the 
form Eb,{T)/Ej = 1 — aln(,/3f),co'c/27r) and (jJr in the expression (31) then becomes replaced by ER{T)/h. For very 
high energy scales, we check that Er{T) would c;onvcrgc to the bare value Ej in the Hamiltonian H. 

These two facts exemplify the asymptotic freedom where microwave light and two-level system almost disentangle. 
Note the parallel between temperature and driving effects in the spin susceptibility close to the confinement frequency. 
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